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ABSTRACT 

The dynamical properties of six SDSS quasar pairs at z ^ 0.8 are investigated. 
The pairs have proper transverse separation R± ^ 500 kpc, and velocity difference 
along the line of sight AV r ^ 500 km/s. If they are bound systems their dynamical 
mass can be evaluated and compared with that of host galaxies. Evidence is found 
of an excess of the former mass with respect to the latter. This suggests that these 
quasar pairs are hosted by galaxies with massive dark halos or that they reside in a 
group/cluster of galaxies. 
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1 INTRODUCTION 

Quasars (QSO) are rare and sho rt-lived objects (e.g., 
iMartinil 12004 iHopkins et ail 120051 ). nevertheless a num- 
ber of associations of QSOs have been discove r ed in 
the last decades (e.g.. IShaverl Il984l : iDiorgovskil Il99ll : 
I Zhdanov fc Surdeil l200ll ). The study of these systems is 
important in the understanding of the evolutionary his- 
tory of galaxies w ith cosmic time and the mechanism of 
QSO ignition (e.g., Di Matteo. Springel. fc Hernquistll2005l ; 



iForeman. Volonteri. fc Dottil 120091 . Particular interest has 
been dedicated to binary QSOs, i.e., two QSOs that reside in 
the same galaxy and that are characterise d by the presence 
of double systems of emis sion lines (e.g., iBoroson fc Lauen 
120091 ; I Rosario et al.l l201ll ) . These systems are thought to 
form in the last s tages of a major merger event (e.g., 
IColpi fc Do tti 2009, and references therein). 

The search of QSO pairs (QSOP) at scales from tens 
to hundreds of kiloparsecs in large surveys was mainly fo- 
cused on_the_Juyestigaton properties 
(e.g., iHennawi et al.l [2010l : IShen et al.ll201ol ) and in partic- 
ular on the excess, with respect to the large scale extrap- 
olation, found at separations of tens o f kilop arsecs (e.g., 
Hcn nawi et al ] l200d : iMvers et all 120071 120081 ). The study 
of the clustering allows us to estimate the bound mass 
of the structures inha bited by QSOs (e.g., ICroom et al.l 
l2005l : IShen et al.ll2010l ). but little attention has been given 
thus far to the study of the dynamical properties of sin- 
gle QSOPs that, if isolated, are dominated by the mass 
of their host galaxies (e.g., iMortlock. Webster, fc Francis! 
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ll999l : lBrotherton et al.ll 19991 ). Although the Cold Dark Mat- 
ter models of galaxy formation predict that QSOs, and in 
particular QSOPs, reside preferentially in particularly rich 
envir onments (e.g.. lEfstathiou fc~R ccs 1988; [Hopkins et "all 
2008), some observational evidence shows that QSOPs could 
be is olated systems (e.g., iFukugita et al.l [20041 ; I Boris et al.l 
120071 ). 

In this paper we look for QSOP s in the Sloan Digital 
Sky Survey (SDSS; I York et al-lBoOol ). with the goal of re- 
constructing the systemic dynamics of the pairs. We found 
six QSOPs at redshift < 0.8, for which the measurement of 
[OIII] lines allows us to pursue this study. In ijUwe describe 
our sample. 33 deals with measurements of radial velocity 
differences. In f|4] we compute virial masses and compare 
them with those of the host galaxies. We investigate the 
QSOP environment in Implications of our results are 
discussed in Sj6] 

Throughout this paper we consider a concordance cos- 
mology with Ho = 70 km/s/Mpc, Q m — 0.3, and Ha = 0.7. 



2 THE QUASAR PAIR SAMPLE 

We investigate the cat alogue of spectro scopically confirmed 
QSOs constructed by ISchneider et al] J2010l ) on the basis 
of the SDSS DR7 (|Abazaiian et all 12009) which contains 
~ 100, 000 objects. We select as pairs two QSOs that have 
proper transverse separation R± ^ 500 kpc, and radial ve- 
locity difference AV r ^ 500 km/s, as based on SDSS red- 
shifts. 14 pairs that satisfy the above criteria are found in 
the redshift range 0.5 < z < 3.3, with luminosities between 
M v ~ -22 and Mv ~ -25. 

Since we are interested in the dynamical properties of 
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these systems we also require that the forbidden [OIII] lines, 
which are used to measure the systemic velocity of the 
QSOs (see ij3|, are present in the SDSS spectra. This im- 
plies that the candidate QSOPs are at redshifts below 0.8. 
With this additional condition we obtain a list of six pairs 
of radi o quiet QSOs (see Ta ble [TJ , five of them considered 
also by iHennawi et al l ^OOd ). 

The probability that they are chance superpositions is 
rather low. In fact, searching for QSOPs in a random sam- 
ple generated with the redshift permu tation method (e.g., 
lOsmerl [l98ll ; IZhdanov fc Surdeil l200lh . which consists of 
maintaining the position of the QSOs fixed, but to randomly 
permute the redshift, we expect to find ~ 0.4 such pairs com- 
pared to the 6 observed. Note that in this new sample most 
of the correlations between objects are destroyed, but the 
angular correlation between QSOs is preserved, so the re- 
sult can be considered as an upper limit for the number of 
chance QSOPs. Thus we assume that all these QSOP are 
physically associated. 

We can exclude that these QSOPs are gravitational 
lens images because: significant differences in the spec- 
tra of the two QSOs are apparent (see Figure [JJ), wide 
separation (AO > 3 arcsec) len sed QSOs are quite rare 
i|Kochanek. Falco, fc Munozlll999h . and there is no evidence 
in SDSS images for luminous galaxies in the foreground of 
QSOPs that could act as a lens. 



3 VELOCITY DIFFERENCES FROM [OIII] 
LINES 

We can also exclude the possibility that the velocity differ- 
ences can be related to the Hubble Flow and therefore mea- 
sure the physical distance of th e pairs. In fact w e verif y that, 
under this hypothesis, in the ISchneider et"aH (|201Ch QSO 
catalogue there are 35 pairs with R± ^ 4 Mpc and physical 
radial separation Rn ^ 4 Mpc. Assuming that the 29 sys- 
tems with R± > 0.5 Mpc are homogeneously distributed, 
we expect ~ 1 with R± ^ 0.5 kpc, while 6 are found. 

It is well known that the redshifts of QSOs de- 
rived from emission lines of various elements can dif- 
fer by as much as 1,000 km/s (e.g ., iTvtler fc Far] 1 19921 ; 
Bon ning. Shields, fc Salviander|[2007f) . Therefore the most 
reliable estimate of the systemic velocity of the QSOs is 
obtained from the measurements of narr ow forbidden lines, 
such as [OIIIJa4949 and [OIII1a50 07 (e.g.. iNelson fc Whittle! 
1 19961 ; lNelsonll2000l : iBorosorJ 12009 ; iHewett fc Wildll2010h . 

We evaluate the baricentres of the lines considering the 
flux above various thresholds with respect to the peak flux 
(see Figure [2}. We take the line position to be the median of 
the individual measurement of the baricentre, and the cor- 
responding uncertainty is given by their interquartile range. 
The redshifts and the radial velocity differences that result 
from these measurements are reported in Table [2] 



4 THE MASS OF QSO PAIRS 

Assuming that the QSOPs form bound systems and thus 
that the velocity difference measured is due to the mutual 
interaction between the two QSOs, we can infer the dynam- 
ical mass through the virial theorem: 
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Figure 1. The SDSS spectra of the two QSOs of each pair 
(background red solid line, foreground blue dotted line). Data 
are smoothed with a 10 A boxcar filter. In order to make visible 
the features of the spectrum of the foreground QSO of QP01, its 
flux is increased by a factor of 10. 



AV 2 R 
G 



(1) 



where AV is the relative velocity of the two components, R 
their separation, and G the gravitational constant. For cir- 
cular orbits, it is possible to calculate the radial component 
of the relative velocity (AV r ) from the redshift difference 
(Az). One has: 

G 



M vi . 



C 



( cAzV 
Vl + J 



(2) 



where c is the speed of light, R± the proper transverse sep- 
aration of the pair (see Table [TJ, and the factor C depends 
only on the inclination angle of the orbital plane t, and on 
the phase angle (f> and is given by: 



C" 



(sin sin i) 2 x \J sin 2 i 



■ + cos 2 m cos 2 b 



(3) 



The average values of C is (C) = 3.4 and the minimum value 

is Grain — 1 - 

In Table[2]we report for each QSOP the minimum virial 
mass (M v ir(min), corresponding to C = 1), which repre- 
sents the minimum mass of the system to be bound. In the 
case of QP06, since there is no significant difference of radial 
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Table 1. Properties of selected QSOPs. z is the redshift from Schneider ct al My is the absolute magnitude in V-band, AO 

and R± are the angular and proper transverse separation, and AV r is the radial velocity difference derived from the redshifts given in 
the catalogue. 







QSO A 


QSO B 








ID 


name 


z 


M v 


z 


M v 


Ad 


R± 


AV r 








[mag] 




[mag] 


[arcsec] 


[kpc] 


[km/s] 


QP01 


SDSS J0117+0020AB 


0.6122 


-22.38 


0.6130 


-24.65 


44 


300 


149 


QP02 


SDSS J0747+4318AB 


0.5010 


-22.76 


0.5012 


-22.61 


9 


56 


40 


QP03 


SDSS J0824+2357AB 


0.5356 


-23.19 


0.5365 


-23.19 


15 


94 


176 


QP04 


SDSS J0845+0711AB 


0.5363 


-23.48 


0.5373 


-23.20 


62 


393 


195 


QP05 


SDSS J0856+5111AB 


0.5425 


-22.81 


0.5434 


-23.59 


22 


139 


175 


QP06 


SDSS J1249+4719AB 


0.4375 


-23.09 


0.4382 


-22.63 


79 


446 


146 



Table 2. Radial velocity difference and virial mass of the QSOPs. z n is the redshift measured from the [OIII] narrow emission lines; 
AV r (n) is the corresponding radial velocity difference; M v i r (min) is the minimum virial mass compatible with the uncertainties of the 
measure of AV r (n). 



ID 


z n (A) 


z„(B) 


AV r (n) 


M v i r (min) 








[km/s] 


[10 12 


M Q ] 


QP01 


0.61142±0.00078 


0.61341±0.00001 


370±171 


2.8- 


-20.4 


QP02 


0.50108±0.00003 


0.50174±0.00001 


132± 7 


0.2- 


- 0.3 


QP03 


0.53527±0.00009 


0.53678±0.00002 


295± 21 


1.6- 


- 2.2 


QP04 


0.53509±0.00015 


0.53754±0.00002 


478± 35 


17.9- 


-24.0 


QP05 


0.54322±0.00003 


0.54239±0.00003 


161± 9 


0.7- 


- 0.9 


QP06 


0.43861±0.00045 


0.43859±0.00001 


4 +94 







velocity, we cannot estimate its virial mass. In this case we 
are probably observing the pair orbit nearly face on. 

It is of interest to compare these M v i r (min) with the ex- 
pected total mass of the pair based on the mass of their host 
galaxies. Acc ording to available mea surements of QSO host 
galaxies fe.g.. lKotilainen et al-llioO^ . and references therein) 
it is found that their mass changes little with redshift. The 
typical range of host mass, based on the galact ic luminosity, 
for o bjects at z < 1 is ~ 0.3 - 1.3 x 10 12 M Q |Decarli et al I 
|2010| . and references therein). 

While for three QSO pairs (QP02, QP03, and QP05) 
their M v i r (min) is consistent with that expected by the 
typical host galaxy masses, in two cases (QP01 and QP04) 
the minimum virial mass is substantially larger than that of 
their host galaxies (see Table [2]). If one assumes the average 
value of C ((C) = 3.4) instead of its minimum, then the 
above cases are further strengthened and also QP03 would 
exhibit a significant mass excess. For the whole (small) sam- 
ple the median value for the M„j r is 6.5 x 1O 12 M0. 

A possible explanation for this mass excess is that QP01 
and QP04 belong to a group or a cluster of galaxies. In this 
case in fact the measured velocity difference depends on the 
overall mass distribution. In the next session we investigate 
this possibility. 



5 QSO PAIRS' ENVIRONMENT 

We searched the SDSS i-band images for an overdensity 
of galaxies that could justify the mass excess discussed 
above. The SDSS magnitude limit in this band is 21.3 mag 
|York et all 120001 1, thus it allo ws us to reach ~ (M* + 1), 
where M* = -20.5 at z = 0.5 jWolf et all 12003), therefore 



Table 3. Environment of QSOPs. n(bkg) is the density of galaxies 
in the region between 2 Mpc and 4 Mpc, N(< 0.5 Mpc) is the 
number of galaxies in the inner 500 kpc, and n(< 0.5 Mpc) is the 
corresponding density. M/L is the minimum mass— to-light ratio 
that could have a galaxy cluster detected on SDSS i-band images 
(3cr over the background). The associated uncertainties represent 
the la statistical fluctuations. 



ID 


n(bkg) 


N(< 0.5 Mpc) 


n(< 0.5 Mpc) 


M/L 




[arcmin -2 ] 




[arcmin -2 ] 




QP01 


1.2±0.2 


5 


1.4±0.6 


> 30 


QP02 


1.4±0.1 


8 


1.4±0.5 


£5 


QP03 


2.1±0.1 


26 


4.7±0.9 




QP04 


1.9±0.1 


15 


2.7±0.7 


^ 100 


QP05 


1.8±0.2 


5 


0.9±0.4 


5 


QP06 


1.9±0.2 


12 


2.3±0.7 





these images permit us to detect only the bright part of the 
galaxy luminosity function. 

The galaxy search was performed using SExtractor 
|f3ertin fc Arnout3[l99 6) on the SDSS images in an area of 
4 Mpc around each pair. The threshold limits for the detec- 
tions is set at 1.5 times over the rms of the background, and 
we classified as galaxy all the sources with the STARCLASS 
parameter lower than 0.2. The number of galaxies in the 
fields (see Table [3} is consistent with the e xpectation from 
the st udy performed up to I = 24 mag by IPostman et al.l 
(1998) on a region of 4° x 4°, and the number of galac- 
tic stars with the p rediction of the TRILEGAL packag43 by 
iGirardi et all l|2005h . 

1 http://stev.oapd.inaf.it/cgi-bin/trilegal_L4 
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Figure 2. The normalised SDSS spectra of the two QSOs of 
the pairs QP03 in the region around the emission line [OIII] at 
A = 5007 A. The central value (red vertical lines) is the median 
value of the baricentre calculated considering different thresholds 
of line peak flux (i.e., 30, 40, 50, 60, and 70%, blue dotted hori- 
zontal lines). 
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Figure 3. Expected central density of galaxies in SDSS i-band 
images as a function of the mass of the clusters for various value 
of the mass— to-light ratio (filled lines). Black square points are 
the QSOPs for which no significant evidences of an overdensity 
of galaxies are present. For these systems we assume that, if a 
cluster is present, it must have a density < 3cr the variation of 
the background. The red square point indicates the pair QP03, for 
which SDSS images highlight a significant overdensity of galaxies 
in the first 500 kpc from the QSOP. 



do not find indications for overdensities larger than 3 times 
the variation of the background, thus, to explain the min- 
imum virial masses of the pairs, these systems require a 
mass-to-light ratio M /L > 5-100 (see Table |3j . Note that 
these values are comparable with those reported in vari- 
ous studies on dyna mical properties of galaxy clusters (e.g., 
IPopesso et alJliooBI . and references therein). 



In order to highlight a possible overdensity around the 
QSOPs, we compute the number of galaxies in annuli of 
500 kpc radius, starting from the centre of each pair. We 
then compare the galaxy density in the first 500 kpc with 
that in the region between 2 Mpc and 4 Mpc, assumed 
as background. These values are reported in Table [3] Only 
QP03 shows a significant overdensity of galaxies. In the other 
cases there is no evidence for a galaxy excess above the back- 
ground by more than 3a. 

We evaluate the expected density of galaxies brighter 
than the SDSS luminosity limits (i ~ 21.3 mag) if a 
cluster of mass M to t = M v i r (min) were associated with 
the QSOPs. We assume that the galaxies of the cluster 
follow th e Schechter lumin osity function with parameters 
given by IWolf et al l (|2003h . and that the galaxies are dis- 
tributed according to a King profile with a virial radius cal- 
cu lated from the virial mass following the relations reported 
bv lGirardi et~ai"1 (|l998l ). 

We compare the expected galaxy density with that ob- 
served in SDSS images (see Figure[3]). In all cases but one,we 



6 CONCLUSIONS 

The analysis of the properties of 6 low redshift QSO pairs 
has shown that in at least two cases the dynamical mass of 
the pair exceeds, by a factor > 10, that expected from their 
host galaxies. A possible explanation of this excess is that 
the QSO host galaxies are surrounded by dark matter ha- 
los w ith masses similar to th ose found in massive ellipticals 
(e.g.. lNapolitano et al.|[2009l ). Alternatively the observed ve- 
locity differences could be due to the presence of a cluster or 
a group of galaxies associated with the QSO pairs. An anal- 
ysis of SDSS i-band images shows evidence for a significant 
overdensity of galaxies in only one case. For the other sys- 
tems a lower limit to the mass-to-light ratio was determined 
at M/L > 5-100 for galaxy clusters with masses equal to 
the virial masses of the pairs. 

In order to strengthen the evidence of a mass excess, 
we can consider a lar ger sample given by th e lists of already 



we can consider a larger sample given by th e lists ot already 
known QSO pairs l|Schneider et al.l l2010l : iHennawi et al.l 
l200dl2010l;lMvers et al.ll200Sl ). Mo"st of these systems are at 
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z > 0.8, excellent instrument capabilities are thus required 
to perform these studies. 
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